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ABSTRACT 

The close active binary HR 5110 was observed at six epochs over 26 days using a global VLBI array at 
15.4 GHz. We used phase-referencing to determine the position of the radio centroid at each epoch with 
an uncertainty significantly smaller than the component separation. After correcting for proper motion and 
parallax, we find that the centroid locations of all six epochs have barycenter separations consistent with an 
emission source located on the KIV secondary, and not in an interaction region between the stars or on the 
F primary. We used a homogeneous power-law gyrosynchrotron emission model to reproduce the observed 
flux densities and fractional circular polarization. The resulting ranges of mean magnetic field strength and 
relativistic electron densities are of order 10 G and 10 5 cm -3 respectively in the source region. 


1. INTRODUCTION 

Late-type stars in short-period binary systems often exhibit 
signposts of enhanced magnetic activity. These include per¬ 
sistent photospheric spots, ultraviolet chromospheric emis¬ 
sion, x-ray em ission from hot coronae, an d non-thermal radio 
emission (e.g. Guinan & Gimenez 1993 ). The conventional 
interpretation is that tidal effects enforce spin-orbit coupling, 
and the consequent rapid rotation drives a robust magnetic dy¬ 
namo on one or both components of the system (e.g. Schrijver 
|& Zwaan|1991| ). 

Close active binaries are categorized by whether one of both 
companions fill their respective Roche lobes. The epony¬ 
mous RS CVn sys tems have_both components detached and 
no mass transfer ( [Hall et aL||1978| , while Algol systems 
are semi-detached witlTone component (typically the cooler, 
more evolved star) filling its Roche lobe, leading to episodic 
mass transfer. Finally, contact binaries have both components 
within each other’s Roche lobes, and continuously transfer 
mass. Curiously, radio luminosities^ do not differ greatly 
among these three classes (Umana et al.| |1998|) , while the 
mean X-ray luminosity of Algol-class systems is 3-4 times 


t y 

lower than RS CVn’s ( [Singh et al.|1996| ). This suggests that 
dynamical processes such as mass transfer and accretion disks 
may be decoupled from emission processes related to en¬ 
hanced magnetic fields. 

An important open question for these systems is whether 
the magnetically-driven enhanced emission arises primarily 
from a single active component or from a magnetic interac¬ 
tion region between the stars. Very-long-baseline radio inter¬ 
ferometry (VLBI) provides a powerful method to address this 
question, since for nearby (d < 100 pc) systems, it can probe 
the structure of the radio coronae at spatial scales smaller than 
the binary orbit. A number of close active binaries have been 
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ever, locating the radio structure within the binary system re¬ 
quires high-precision astrometric techniques such as rapidly 
switching between the target star and one or more extragalac- 
tic phase reference sources with precisely determined posi¬ 
tions. Multi-epoch observations over several years are needed 


in order to simultaneously solve for the astrometric param¬ 
eters (parallax, proper motion, orbital motion) with sub-mas 
accuracy. 


The fir_st measurements of this type were made by [Lestrade 
et al. (1993), who showed that the radio centroid of Algol ^s 


inner binary moved with the KI V secondary. Subsequent as¬ 
trometric VLB I studies of Algol ( [Peterson et al.|2010| [Peter¬ 
son et al.||2011| ) confirmed the K star association. They also 
showed that during flares, the radio structure consisted of a 
large coronal loop approximately one stellar diameter height, 
with a base straddling the sub giant and oriented toward the B 
star. Other recent astrometric VLBI studiesof close active bi¬ 
naries include UX Arietis (Peterson et al. 2011) , and IM Peg 
( [Ransom et al.|2012j ). In all three cases, the radio emission is 
associated with the more evolved star, a K subgiant or giant. 

In this paper we report on an astrometric VLBI study of a 
fourth close active binary. HR 5110 (BH CVn, HD 118216) 
is a short period (2.61 day) binary composed of a F2 IV pri¬ 
mary and K2 IV secondary. The or bit is circular and is ori - 
ented nearly face-on to the observer ( [Eker & Doherty |l 987 1 . 
This orientation is well-suited to investigate whether the ra- 
dio emission is associated with a single component or with 
an interaction region. It also provides a convenient test of ac 
five region association with polar spots (Huenemoerder et al. 
2009[ ), assuming the magnetic and rotational axes are parallel 


and normal to the orbital plane, and hence oriented toward the 
observer. 

Although no spots have been directly detected on HR5110, 
Little-Marenin et al. (1986) found that the IR spectrum can 


be best-fit by assuming a persistent spot on the K secondary 
with a filling factor ~0.25. J.t has been categorized both as 
an RS CVn syst em ([Hall et ak||1978[ ) and an Algol system 
(Little-Marenin et al. 1986 ), although UV observations have 
foundno direct evidence of a accretion flow or disk. However, 
Little-Marenin et al. (1986 ) argue that the secondary must fill 
its Roche lobe based on the spectral classification and its close 
binary separation. 

There have been several estimates of the size and location 
of the active r egion in HR5110. Previous VLBI studi es of the 
radio corona |Mutel et al.|1985[ |Ransom et al.|2003] ) found a 
‘core-halo’ structure, with a compact core size of order the K- 
star radius and an more extended halo somewhat larger than 
the orbital separation. This is similar to the results of Umana 
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et al. ( 1993| ), who measured the radio spectrum of HR5110 
betweeiTl.4 GHz and 15 GHz during several epochs of mod¬ 
erate activity. They found that the spectral energy distribution 
could be fit with a two component gyro-synchrotron model 
consisting of higher density core with a size about one-half 
the K-star radius and a lower density extended component ap¬ 
proximately the size of the binary system. 

Observations at other wavelengths provide additional con¬ 
straints on the active region’s size and location. Little 
|Marenin et al.| ( |1986] ) found that UV emission lines are best-fit 
with the K secondary’s chromosphere and transition region, 
hence locating UV flares on the K star. |Mullan et al.| ( |2006| ) 
modeled UV flare observations with a compact "active loop 
whose length was smaller than 0.42 stellar radii. However, 
|Graffagnino et al. | ( | 19951 ) argued that the active region may be 
much larger. They fit ROSAT X-ray data from a large out¬ 
burst lasting more than three days. They found that a very 
large emission region, on the scale of the binary separation, 
is required to fit the observed light curve. They suggest that 
the flare may have been driven by magnetic reconnection in 
the interbinary region between the stars. Assuming the ra¬ 
dio, UV, and X-ray emission regions share a common active- 
region site, an open question is: Where in the binary system 
does this region lie? 

We observed HR5110 at multiple epochs using an intercon¬ 
tinental phase-referenced VLB I array. The observations were 
made at a frequency of 15 GHz, which provided twice the 
angular resolution of any previous VLB I study of this sys¬ 
tem. The primary goal of the observations was to determine 
the location of the radio emission within the binary system. 
We also modeled the observed emission properties (flux den¬ 
sity, circular polarization, angular size) using a power-law gy- 
rosynchrotron emission model to derive constraints on plasma 
parameters in the source region. 


2. OBSERVATIONS AND DATA ANALYSIS 

HR5110 was observed at six epochs of ten hours each with 
a global VLB I array consisting of the ten 25 m telescopes of 
the Very Long Baseline Array (VLBA^]and the 100 m Effels- 
berg telescope near Bonn, Germany. We scheduled observa¬ 
tions with a cadence that allowed sampling at approximately 
evenly spaced orbital phases. The receivers sampled 64x2 
MHz channels using 2-bit sampling centered at 15.4 GHz in 
dual polarization mode. The synthesized beam was 0.4 mas 
in its smallest dimension which is approximately one fifth the 
binary’s major axis (Table [TJ. 

All epochs were calibrated and imaged using the NRAO 
Astronomical Image Processing Software package (AIPS, 
Grei sen|2003| . Standard VLBI amplitude and delay-rate cor¬ 
rections were applied, followed by the transfer of complex 
gain corrections from the primary calibrator (J1308+3546) 
to the target (HR5110) and secondary phase calibrator 
(J1310+3220) (Diamond 1995j ). We assumed a flux density 
of 0.6 Jy for J1308+3546. After calibration, we imaged all 
sources using Fourier inversion and numerical deconvolution 
(AIPS task IMAGR using the Clark CLEAN algorithm). We 
used the nodding phase-referencing scheme (Lestrade et al. 
T990|, cycling at three minute intervals between J1308+3546 
and HR5110. We observed J1310+3220 once every 15 min¬ 
utes to check the stability of the phase referencing scheme. 


1 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 


The nodding phase-referencing technique requires accurate 
interpolation of phases from calibrator scans to target scans. 
This proved to be problematic, both because of the the rela¬ 
tively large angular separatior0 between phase calibrator and 
HR5110 (5.5°), but also because all epochs were scheduled 
during summer months, which often results in large short- 
timescale phase fluctuat ions caused from wet-component tro¬ 
pospheric variations ( Pradel et al.||2006 ). Although we able 
to successfully transfer phase solutions from calibrator to tar¬ 
get scans, it required careful inspection and flagging of exces¬ 
sive phase jumps (AIPS task SNFLG). We flagged all time- 
baseline visibilities for which the phase difference between 
adjacent 30 second calibrator phase solutions exceeded 60°, 
indicating a possible phase wrap ambiguity over the three- 
minute nodding cycle. This resulted in more than half the 
visibilities in each ten-hour observing epoch being flagged. 

In spite of this significant data loss, the positional accuracy 
of the radio centroid positions was only modestly degraded. 
To establish the positional accuracy at each epoch, we mea¬ 
sured variations in the measured position of the secondary 
phase calibrator (J1310+3220, angular separation 3.3° from 
J1308+3546) compared with its ICRF position. At all epochs 
but epoch A (2012.495), the measured position was within 
±0.1 mas of its ICRF position. This is c onsis tent with the 
extrapolated uncertainty estimate of Prad el et al.| ( |20Q6] ), who 
studied the astrometric accuracy of VLB A phase-referenced 
observations as a function of source-calibrator angular sepa¬ 
ration and declination for a range of tropospheric conditions. 

The expected astrometric accuracy of phase-referencing 
can be estimated by linearly extrapolating the RMS position 
uncertainty trend given in Figure 2 of |Pradel et al.| |2006| . 
For a source at declination 8 = ±25°, and source-phase refer¬ 
ence angular separation A0 degrees, at an observing frequency 
15 GHz and with a significant wet tropospheric component 
this can be written as 


o(A0) ~ 30 inns x A0 

where a is the positional uncertainty in micro-arcsecond for 
VLB A observations. Hence, the expected total RMS posi¬ 
tional uncertainty at 3.3°separation (primary-secondary cali¬ 
brator) is 0.1 mas, in excellent agreement with our observa¬ 
tions, and 0.17 mas at 5.5°separation (primary-target sepa¬ 
ration). We have used this estimate for the formal position 
uncertainties at epochs B-F. For epoch A, the peak-to-peak 
position jitter of the secondary calibrator was 0.2 mas, possi¬ 
bly resulting from an excessively wet troposphere at several 
telescopes. Scaling this uncertainty to 5.5°separation resulted 
in an estimated uncertainty of 0.33 mas for the position of 
HR5110 at Epoch A. 

The radio centroid could move by up to 0.7 mas over the 
course of each ten hour observing epoch, if the source were 
located at the KIV secondary. This position shift would sig¬ 
nificantly smear out the source brightness in maps made us¬ 
ing the entire 10-hour observation. The motion can be cor¬ 
rected by adjusting the map phase center (using AIPS task 
CLCOR) with a model of the source location within the orbit 
versus time. This is problematic since the source position is 

2 Previous phase-referenc ed VLBI observations of HR5110 jLes trade et al.| 
|1999[ |Ransom et al.||2003) used several angularly closer phase calibrator 
sources, but at lower frequencies (5, 8 GHz). To check their suitability at 
15 GHz, in April 2012 we made ’snapshot’ VLB A observations of several an¬ 
gularly nearby compact sources. We found that the previously-used sources 
were heavily resolved at 15 GHz and therefore unusable, and that the closest 
suitable compact source was J1308+3546. 
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TABLE 1 

HR5 110 PROPERTIES AND ORBITAL ELEMENTS 


Parameter Symbol 


Value 


Ref. a 


Fixed values 


Parallax 

n 

22.21 mas 

1 

Spectral type 


F2IV + K2IV 

2 

Primary mass 

m A 

1.5 M q 

2 

Secondary mass 

m B 

0.8 Mq 

2 

Primary radius 

r A 

2.6 Rq (0.27 mas) 

2 

Secondary radius 

r B 

3.4 R q (0.35 mas) 

2 

Reference epoch [phase 0.0] b 

T 

2445766.655 

2 

Eccentricity 

e 

0.00 

2 

Inclination 6 

i 

171.1° 

2 

Semi-major axes 

a\ ,«2 

0.017 AU (0.38mas), 0.032 AU (0.71 mas) 

2 

Component separation 

a\ +<3 2 

0.049 AU (1.09 mas) 

2 

Period 

P 

2.613214 day 

3 

Derived values 

Longitude of ascending node 6 

a 

89° ± 10° 


R.A. proper motion 

Pa COS 5 

85.62 ± 0.07 mas yr -1 


Dec. proper motion 

P6 

-9.68 ± 0.07 mas yr -1 


Fiducial d R.A. 

OCo 

13:34:47.75949 + 0.00008 


Fiducial d Declination 

5o 

37:10:56.7605 + 0.0006 



a References: (1) Lestrade et al. (1999), (2) Eker & Doherty (1987), (3) Mayor & Mazeh| fT~987) 
b Secondary is in conjunction (in front) (Eker & Doherty|1987|i 
c Assumes radio centroid is centered on KTsubgiant (see text). 
d Center of mass position at fiducia l epoch (JD = 2449028.000) 

e Inclination from |Eker & Dohertylfl987) , but corrected for CW rotation (c.f section 3). 


not known a priori! In our case, this correction was unneces¬ 
sary, since the flagging described above resulted in using only 
^ 3 hr contiguous unflaggged time intervals at each epoch. 
These were typically near the end of each 10-hour observa¬ 
tion. Because the source motion on this timescale is less than 
a resolving beam, there was no significant smearing effect. 

In addition to the VLB I observations in this paper, we also 
calibrated and imaged five unpublished phase-referenced ob¬ 
servations of HR5110 from the public VLB A archive. Al¬ 
though these observations lack the spatial resolution and tem¬ 
poral sampling cadence to determine the position of the radio 
source within the binary orbit, the measured positions, along 
with the published position of Lestrade et al. ( 1999 ), allowed 
us to make a more precise determinationof the proper motion 
of HR5110. 


3. RESULTS 

We determined new proper motions for HR5110 (Table [T} 
by combining our 6-epoch position determinations with the 
previously published position of |Lestrade et al. ( |1999| ) and po¬ 
sitions that we determined from archival VLB A observations, 
as listed in Table [2] Our calculated right ascension proper 
motion is in good agreement with the Lestra de et al.| (1999) 
value, while the declination proper motion isThree"standard 
deviations from their value. 

Radio contour maps of HR5110 at all six epochs are shown 
in Figure [I] Previous VLBI observatio ns of HR5110_found a 
core-ha lo source morpholog y at 5 GHz (Mutel et al. 1985 ) and 
8 GHz ( [Ransom et al.|2003| ). The 15 GHz maps are somewhat 
amorphous and do not manifest a clear core-halo structure 
even if the long baselines to the more sensitive Effelsberg tele¬ 
scope are excluded. However, we can characterize the spatial 
scale sizes and fluxes by fitting the images with a one compo¬ 
nent Gaussian. Fluxes vary significantly and are likely under¬ 
estimates of the true flux due to the rapidly varying wet tropo¬ 
sphere introducing significant phase errors (Marti-Vidal et al. 


[2010| . All epochs were fit using a single component Gaus¬ 
sian and the resulting Gaussian model parameters are given in 
Table [3] 

We used our mea sured proper motions, along with the par¬ 
allax determined by Lestra de~et al.| (|T999), to convert the ob¬ 
served radio centroids to displacements relative to the binary 
center of mass. These displacements were used to solve for 
the co-rotating location of the radio source with respect to the 
binary’s center of mass. To do this, we assumed that the ra¬ 
dio emission at all epochs is located at a fixed point in the 
co-rotating frame. We solved for this location by finding the 
point in the co-rotating frame that minimizes the sum of the 
squared distances (weighted by position uncertainty) from the 
center of mass to the radio centroids. Figure l] shows the % 2 
minimization surface as a function of the radial displacement 
from the center of mass and an azimuthal angle measured 
north through east in the sky plane. Table[4]lists the positional 
offset at each epoch between the measurea radio centroid and 
the best-fit fixed point in the co-rotating frame. 

The best fit radial displacement is nearly the same as the 
semi-major axis of the K2 IV secondary, rather than the F2 
IV primary or an interaction region midway between the stars. 
Since the binary’s orientation on the sky is not known, the 
radio centroid azimuthal angles cannot be directly compared 
with the location of the component stars. However, given the 
close agreement between the model’s radial displacement and 
the K star’s semi-major axis, we make the plausible assump¬ 
tion that the radio centroids are located on or near the K star 
rather than a random co-moving location on the K star’s or¬ 
bit. Figure [3] shows the radio centroids and star positions in 
co-moving sky plane coordinates, while Figure |4]is in the co¬ 
rotating frame, assuming the radio centroid-K star identifica¬ 
tion described above. 

The identification of the radio centroid with the K IV sec¬ 
ondary allows us to determine new values for two orbital el- 
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TABLE 2 

J2000 Heliocentric Positions 


JD 

Year 

Rt. Asc. (13:34:+) 

Decl. (37:10:+) 

Phase Reference. 3 

Ref. b 


Archival and previously published 

2449028.000 

1993.108 

47.75947 ± 0.00003 

56.7601 ± 0.00070 

OP326 

1 

2451270.833 

1999.250 

47.80353 ± 0.00020 

56.6997 ± 0.00200 

J1340+3754 

2 

2451276.833 

1999.266 

47.80374 ± 0.00020 

56.7003 ± 0.00200 

J1340+3754 

2 

2451631.831 

2000.239 

47.81060 ±0.00004 

56.6915 ± 0.00060 

J1317+3425 

3 

2454560.812 

2008.258 

47.86794 ± 0.00003 

56.6129 ± 0.00030 

J1340+3754 

4 

2454562.813 

2008.264 

47.86799 ± 0.00003 

56.6147 ± 0.00020 

J1340+3754 

4 

This paper 

2456108.458 

2012.495 

47.89844 ± 0.00002 

56.5721 ± 0.00034 

J1308+3546 

5 

2456118.438 

2012.523 

47.89852 ± 0.00003 

56.5722 ± 0.00017 

J1308+3546 

5 

2456119.427 

2012.525 

47.89865 ± 0.00001 

56.5725 ± 0.00017 

J1308+3546 

5 

2456122.417 

2012.534 

47.89867 ± 0.00001 

56.5726 ± 0.00017 

J1308+3546 

5 

2456133.385 

2012.564 

47.89883 ± 0.00001 

56.5728 ± 0.00017 

J1308+3546 

5 

2456134.396 

2012.566 

47.89888 ± 0.00001 

56.5715 ±0.00017 

J1308+3546 

5 


a Positions in reference to the calibrators: OP326 (RA=13:17:36.49418 Dec=34:25:15.9326), 
J1340+3754 (RA=13:40:22.9518 Dec=37:54:43.835), J1317+3425 (RA=13:17:36.4942 

Dec=34:25:15.9 33), J1308+3546 (RA = 13:08:23.70911 Dec = 35:46:37.164) 

b References: 1. Lestrade^t al. (19991, 2. VLBA archive code BG087, 3. VLBA archive code BB120, 
4. VLBA archive S9644, 5. This paper. 


TABLE 3 

Observed Fluxes and Gaussian model parameters 


Epoch 

Code 

<p a 

Stokes I (mJy) 

Stokes V/I 

Diameter (mas) b 

log(7h) 

2012.495 

A 

0.54 

2.6 ± 0.6 

<0.28 

<1.04 

>7.1 

2012.523 

B 

0.37 

7.7 ± 0.8 

<0.08 

0.76 ± 0.20 

7.8 ±0.1 

2012.525 

C 

0.73 

76.5 ± 5.0 

+0.02 ± 0.01 

0.86 ± 0.05 

8.7 ±0.1 

2012.534 

D 

0.88 

91.6 ±0.4 

+0.04 ± 0.01 

0.36 ± 0.05 

9.6 ±0.1 

2012.564 

E 

0.11 

12.2 ±0.9 

+0.18 ±0.01 

0.94 ± 0.20 

7.9 ±0.1 

2012.566 

F 

0.47 

10.3 ±1.6 

<0.12 

1.08 ±0.20 

7.7 ±0.1 


a Orbital phas e computed using the period in Table [I] and adopting phase 0 at secondary con¬ 
junction (Eker & Doherty|1987) . 

b Sizes are the full-width at half maximum geometric means of elliptical single component 
Gaussian fits. 


TABLE 4 

Offsets from best-fit model position 


Epoch 

Aa (a) a 

A5 (a) a 

1993.108 

0.0 (0.0) 

-0.9 (1.3) 

1999.250 

1.1 (0.6) 

-1.0 (0.5) 

1999.266 

1.4 (0.7) 

0.1 (0.1) 

2000.239 

0.7 (1.4) 

0.8 (1.4) 

2008.258 

-0.3 (1.0) 

-0.7 (2.5) 

2008.264 

-0.5 (1.7) 

0.3 (1.1) 

2012.495 

0.6 (3.9) 

-0.1 (0.6) 

2012.523 

-0.2 (0.6) 

0.1 (0.3) 

2012.525 

0.2 (1.5) 

0.0 (0.0) 

2012.534 

-0.1 (0.4) 

-0.4 (2.6) 

2012.564 

0.0 (0.1) 

0.1 (0.8) 

2012.566 

0.0 (0.2) 

0.0 (0.1) 


a Observed - model in mas (std. dev.) 


ements of the binary system. First, since the observed radio 
centroid motion in the center-of-mass system is clockwise on 
the sky (Figure [3j black arrow on secondary orbit), the incli¬ 
nation angle is in the second quadran t (|Heintz|1978|). Hence i 
= 171.1° rather than 8.9° as given by Eker & Doherty](1987). 


Second, the longitude of the ascending node (labeled Q in 
Figure [3]), which was previously unknown, was determined 
by using the epoch of secondary conjunction (Eker & Doherty 


1987] phase 0) to determine the interpolated radio position at 
phase 0.75, i.e. the time of transit through the ascending node. 
We find Q = 89° ±10° where the uncertainty corresponds to 
the scatter in co-rotating radio centroid positions (Figure [4]) 
translated to orbital phase. 

4. DISCUSSION 


4.1 . Active Region Location 

The principal astrometric result of this study is that the radio 
centroids of HR 5110 are located within one stellar radius of 
the co-rotating KIV secondary center at all epochs. The iden¬ 
tification of the radio emission with the secondary component 


such systems viz., Algol (Lestrade et al. 1993, 

Peterson et al. 

2010), UX Ari ([Peterson et al. 2011), and HR! 

3703 (Ransom 


et al. 


20121. Similarly, |Little-Marenin et al~ ( [19861 ) found that 
UV emission lines associated with an active chromosphere 
and trans ition region are located on the KIV secondary. By 
contrast,[Graffagnino et al.|(jl995J) modeled a very large X-ray 
flare from HR5110 and found that the source size was com¬ 
parable with the binary separation. They suggested that the 
interaction region between the stars could be the origin of the 
emission, e.g., t he interbinary magnetic r egion posited for the 
binary CZ CYn (Strassmeier et al. |2011| ). In addition, during 
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Fig. 1.— Radio contour maps of HR5110 at 15 GHz overlaid on the binary system in the center of mass frame of reference. The contour intervals are: 50, 
60, 70, 80, 90 percent of the maximum radio intensity at each epoch. The axes label ticks are in mas relative to the binary center of mass. The superposed 
star positions (green circle= FIV primary, red circle = KIV secondary) are consistent with a least-squares minimization solution for a fixed location of the radio 
centroids (at all six epochs) in the co-rotating frame of the binary system (see section 3). Since the resulting center of mass radial displacement was close to the 
K-star displacement, we identify the radio source location with the K secondary. The grey ellipse is the restoring beam at each epoch. 
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Fig. 2.— Confidence intervals of best-fit solution (center of mass displace¬ 
ment, ascending node angle) for stationary location of radio emission in the 
co-rotating frame. The dotted lines indicate the radial offsets from the center 
of mass for the K2 IV star (red), F2 IV star (blue) and the midpoint between 
the two (green). Contours represent confidence intervals centered on the best- 
fit position (r = 0.70 mas, Q. = 89°). 



Fig. 3.— HR5110 radio centroid locations at six epochs in the co-moving 
center of mass system. Each epoch is color coded: A (yellow), B (orange), C 
(blue), D (cyan), E (magenta), F (black). The large black circle is the orbit of 
the cooler KIV secondary star. Smaller colored circles show the secondary’s 
size and location at the time of the corresponding observation. The dotted 
lines indicate the distance between the radio centroid and the KIV center at 
each epoch. 


Fig. 4.— Radio centroid positions overlaid on the binary system in the co¬ 
rotating frame. The orientation of the stars is based on the close association of 
the best fit radial offset with the K2 IV semi-major axis. Each epoch is color 
coded: A (yellow), B (orange), C (blue), D (cyan), E (magenta), F (black). 
The two circles are the K2 IV star (red) and the F2 IV star (green). 


an extreme radio flare Algol’s emission was c learly displaced 
towards the hotter component (Peterson et al. |2010| ). 

This apparent disagreement concerning the active region 
location poses an interesting conundrum: Are these systems 
driven by magnetic activity on the cooler star or between the 
stars? It is possible that the size and location of the active re¬ 
gion depends on the intensity and duration of the flare. For all 
our VLBI epochs, the observed radio flux density of HR5110 
was significantly les s than the highest recorded flux outburst 


_ignest rec 

0 S > 400 mJy, Feldman & Viner [l980| ). However, the large 
HR5110 X-ray event and the large radio loop observed on Al¬ 
gol represent extraordinarily large outbursts. It is possible that 
moderate flares are located on the cooler secondary but that 
emission regions responsible for the largest flares can extend 
into, or arise from an interbinary region. 

Previously the radio emission was thought to be tied to 
a magnetically active polar spot ( [Elias et al.||1995| ), but the 
VLBI determined radio positions are not tightly grouped on a 
scale smaller than a stellar radius. Polar spots in bin ary sys¬ 
t ems are thought to be long lived, up to several years ( Strass 


meier 2009). By contrast the scatter of radio positions seen in 
Figure [4] occurs over a time period only slightly more than a 
month. In HR5110 the orbital plane is inclined nearly perpen¬ 
dicular to the line of sight. If the magnetic pole of the cooler 
star is roughly aligned with its rotation axis then a putative 
polar spot would be nearly stationary from the observer’s po¬ 
sition. The observed scatter in the HR5110 radio centroids is 
inconsistent with emission from a stationary location associ¬ 
ated with a polar spot This result for HR5110 is very similar 
to HR8703 ( [Ransom et al.|2012| ) where the radio centroids at 
multiple, well-separated epochs are scattered over the entire 
stellar disk. 


4.2. Physical Properties of Emission Region: 
Gyro synchrotron Model 
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Fig. 5 — Comparison between observed and model 15 GHz flux density (colored surface) and fractional circular polarization (CP, black contours) for HR5110. 
Model fluxes were calculated using a uniform box of size 5.8 R & (0.6 mas) at 45 pc emitting gyrosynchrotron radiation for two different magnetic field orientations 
and electron power-law indices: (a) 0 = 60°, 5 = 1.5; (b) 0 = 80°, 5 = 1.5; (c) 0 = 60°, 5 = 2.0; (d) 0 = 80°, 5 = 2.0. Also shown are rectangular boxes delineating 
observed flux densities and fraction CP (with uncertainties) at epochs 2012.525 (C, red) 2012.534 (D, orange), as well as 15 GHz fluxes and CP reported by 
jUmana et al. 1(1993} at epochs 1989.134 (Ul, green) and 1989.178 (U2, cyan). 


To characterize the plasma conditions responsible for the 
observed radio emission, we constructed a simplified coro¬ 
nal model consisting of constant density, mildly-relativistic 
power-law electrons in a uniform magnetic field emitting gy¬ 
rosynchrotron emission. The emitting volume is a cube whose 
size is the average observed VLBI angular size (0.6 mas = 
5.7 Rq). This is clearly an oversimplified model, but allows 
for order-of-magnitude estimates of the electron density and 
magnetic field strength without making unsubstantiated as¬ 
sumptions concerning the magnetic field topology or ener¬ 
getic electron population parameters. 

The model solves for the emergent Stokes I and V flux 
density at the distance to HR5110 (45 pc.). We used the 
expressions of [Robinson & Melrose|(|1984|) for the emission 
and absorption coefficients of gyro synchrotron emission from 
mildly relativistic electrons with a power-law electron energy 
distribution given by, 


where K = No (8 — 1 )£o 5-1 and is 8 the energy index, and with 
a low-energy cutoff Eq - 10 keV. The emergent flux density 
was calculated by integrating the radiative transfer equation, 
summing the contributions along uniformly spaced lines of 
sight through the cube We varied the strength and orientation 
of the magnetic field as well as the density and energy index 
of the relativistic electron population to best-fit the measured 
15 GHz integrated flux density and fractional circular polar¬ 
ization at each epoch. 

Figure [5] shows a surface plot of Stokes I flux density as 
a function of electron density and magnetic field strength for 
four pairs of energy index and magnetic field orientation angle 
0 measured with respect to the observer’s line of sight. Multi¬ 
frequency o bser vations of HR5110 were obtained by Umana 


N e (E) = KE 


-5 


|et al.| ( [T993] ) and | White et al.| ( l990| ). Both reported relatively 
flat spectral indices (a < -0.6) at 15 GHz. For a synchrotron 
process, the spectral index is simply related to the power-law 
electron energy index 8 = 1 — 2a < 2.2, so we chose repre¬ 
sentative values 1.5 and 2.0. The magnetic field orientation is 
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unconstrained by prior observations, but values smaller then 
60° resulted in model circular polarization fractions that ex¬ 
ceeded our observed values, so we show only larger angles. 

For a given Stokes I flux density the solution is highly de¬ 
generate, with constant flux density values occurring along 
contours of nearly constant slope in log(Vg)-log(5) coordi¬ 
nates. These lines have a slope that depends on primarily on 
the energy index, and can be approximated by 


log (Me) 

log(B) 


-0.85 


However, a measurement of circular polarization (CP) can 
break this degeneracy, as shown by the black contours, which 
are labeled by fractional CP. At epochs C and D, the signal-to- 
noise ratio was sufficient to measure significant non-zero cir¬ 
cularly polarized emission (Table [3j, so we used these epochs 
to constrain the model. Observed flux densities and frac¬ 
tional circular polarization are shown as overlaid rectangles 
(red: epoch C, orange: epoch D), where the widths represent 
the uncertainty in each value. For comparison, we also show 
15 GHz flux densities and CP reported by Umana et al. ( |1993| 
at epochs 1989.134 (green, label Ul) and 1989.178 (cyan, la¬ 
bel U2). 

Comparing the model with observed fluxes and CP, we find 
that although there is still a significant parameter degeneracy, 
the present observations favor electron densities in the range 
n e ~ 10 5 cm -3 and magnetic field strengths B~ 10 Gauss, 
but with a large allowable range ( appro ximately db 1 dex). 
By comparison, the Umana et al. ( |1993| observations favor 
somewhat lower electron densities and higher magnetic field 
strengths. Since the source is highly variable, these differ¬ 
ences may reflect different physical conditions in the emission 
region at different epochs. 


5. SUMMARY 

We observed the close active binary HR5110 with a global 
VLBI array at 15.4 GHz at six epochs distributed uniformly 
over the binary’s orbit phase. The primary goal of these ob¬ 
servations was to determine the location of the radio source 


within the binary system. To do this, we used a phase- 
referencing nodding scheme, rapidly switching from the tar¬ 
get star to an angularly nearby extragalactic radio source. 

After correcting for parallax and proper motion, we de¬ 
termined the locations of the radio centroids within the co¬ 
rotating binary system with sub-millarcsecond accuracy. The 
positions were scattered around a fixed location [within mea¬ 
surement uncertainty] in the co-rotating frame of the binary. 
The radial offsets from the center of mass were nearly the 
same as the semi-major axis of the KIV secondary, strongly 
suggesting the KIV secondary as the source of the radio emis¬ 
sion. This identification of the radio source with the cooler, 
more evolved star has been previous reported for sev eral other 
close binary systems , including Algol ([Lestrade et al.| ] 1999; 
|Peterson et al.|201Q|), UX Ari ( [Peterson et al.|2011| ), and IM 
Peg ([Ransom et al.|2012|). 

The clockwise radio centroid motion indicates an orbital 
inclination angl e of 171.1 ° rather than 8.9° as given by |Eker &| 
|Doherty[ ( |l987 ). Additionally, the longitude of the ascending 
node (£2 = 89°) was determined using the identification of the 
radio centroids with the K IV secondary and the previously 


determined time of secondary conjunction (Eker & Doherty 
[1987] ). 


We used the observed properties of the radio emission 
(flux density, angular size, circular polarization) to constrain 
plasma parameters in the source region. The model used a 
simplified geometry consisting of power-law relativistic elec¬ 
trons in a uniform magnetic field radiating gyro-synchrotron 
emission. The derived values of magnetic field strength and 
relativistic electron density are degenerate, but the degeneracy 
is broken when there is a measured fractional circular polar¬ 
ization. For two such epochs, we find a mean magnetic field 
strength of order 10 Gauss and electron densities that vary 
from 10 4 cm -3 — 10 6 cm -3 depending on the electron power- 
law index and magnetic field orientation. 
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